The entropy of an acoustic black hole in Bose-Einstein condensates 
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We compute the entanglement entropy associated to the Hawking emission of a (1 + l)-dimensional 
acoustic black hole in a Bose-Einstein condensate. We use the brick wall model proposed by 't 
Hooft, adapted to the momentum space, in order to tackle the case when high frequency dispersion 
is taken in account. As expected, we find that in the hydrodynamic limit the entropy only depends 
on the size of the box in the near-horizon region, as for gravitational (f + f )-dimensional black holes. 
When dispersion effects are considered, we find a correction that depends on the square of the size of 
the near-horizon region measured in units of healing length, very similar to the universal correction 
to the entropy found in the case of spin-f /2 Heisenberg XX chains. 



I. INTRODUCTION 



When Hawking realized that a Schwarzschild black hole 
emits radiation like a black body with a temperature de- 
termined by its mass [![, investigations focused on the 
connection between the Hawking-Bekestein formula Q 
for the entropy S = jM^A, where A is the area of the 
horizon, and some consistent microscopic counting of de- 
grees of freedom. Generally speaking, there is a large con- 
sensus on defining the entropy by using the von Neumann 
formula S = — Tr(plnp) where p is a density matrix. In 
particular, one can associate the density matrix to the 
sub-state formed by the outside region of the black hole. 
In this case, the entropy measures the degree of entangle- 
ment between the modes in the two sides of the horizon 
Q. Alternatively, entropy can be defined through the 
statistical mechanics of a system in the vicinity of the 
horizon 0, [f| . Remarkably, the two characterizations co- 
incide and agree with the Bekenstein-Hawking formula 
up to the factor 1/4, whose origin is still unknown. It 
should also be mentioned that these microscopic realiza- 
tions of the entropy suffer from ultraviolet divergences, 
that can be cured by introducing a cutoff at around the 
Planck scale, see Q for recent developments and Q for 
a review. 

The lack of experimental evidence for Hawking ra- 
diation is mainly due to the smallness of h and c _1 . 
However, as first noticed by W. Unruh in 1981, in con- 
densed matter physics there are systems that closely 
mimic curved spacetime configurations, and where the 
speed of light is effectively replaced by the speed of sound 
waves, so the suppression of quantum effects can be lifted 
by several order of magnitudes @. In particular, an ir- 
rotational fluid flowing through a device able to accel- 
erate it to supersonic speed can generate a thermal flux 
of phonons that shows the same characteristics of the 
Hawking radiation emitted by a black holes, see [l(| and 
references therein. This possibility was studied in the 



'Electronic address: rinaldi.max@gmail.com 



context of Bose-Einstein condensates (BEC) and many 
other systems, see e.g. [ll|. Although no analog formu- 
lae to the Hawking-Bckcnstein one are known for these 
dumb holes, we expect that entropy can be associated to 
the phonons created via the Hawking mechanism. In fact, 
the acoustic horizon acts as a partitioning screen, which is 
a sufficient condition to create entanglement and, there- 
fore, entanglement entropy. 

In this paper, we would like to address the calcula- 
tion of the entanglement entropy associated to the analog 
of the Hawking radiation created in (l+l)-dimensional 
acoustic black holes in dilute BEC gas. In order to avoid 
typical infrared divergences occurring in this kind of bi- 
partite systems, we confine the region in which we com- 
pute the entropy into a box of size L and located at 
an arbitrary distance e near the horizon. The (1+1)- 
dimensional acoustic black holes was intensively stud- 
ied both analytically [l^ - [l5| and numerically [l6| as it 
might be experimentally realizable. In the limit where 
the wavelength of the modes are much larger than the 
healing length of the gas, one can neglect the high fre- 
quency dispersion typical of this system (the so-called 
hydrodynamic limit). In this case, we expect that the 
entropy is proportional, at the leading term, to ln(L/e). 
This is due to the fact that the mode equation in (1 + 1) 
dimensions is nearly conformally invariant, exactly like 
in the case of a (1 + l)-dimensional gravitational black 
hole [25| . Therefore, the entropy is purely "geometric" 
and arbitrary, in the sense that it cannot depend on the 
parameters of the black hole. In this paper, we verify 
this by employing the brick wall model proposed by 't 
Hooft in see also [l?}- We stress that the entangle- 
ment entropy computed here refers only to the phononic 
radiation produced by the Hawking mechanism, and has 
nothing to do with the thermodynamic entropy of the 
Bose-Einstein condensate, which vanishes. 

The main results of our work concern however the case 
when dispersion is taken in account, and conformal in- 
variance is broken. In fact, the dispersion typically in- 
troduces a high-order differential operator in the mode 
equations, with a prefactor that depends on the healing 
length, in analogy with certain gravitational models en- 
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dowed with modified dispersion relations, see e. g. [18j- 
Because of this term, the entropy is no longer arbitrary 
and it is reasonable to expect that it depends on the heal- 
ing length. Indeed, we find that this is the case, by us- 
ing the brick wall technique introduced by 't Hooft and 
adapted to the momentum space. As far as we know, 
this is the first time that such a method is employed in 
the context of acoustic black holes, and we find that our 
adaptation to momentum space becomes a powerful tool 
to compute the entanglement entropy, especially when 
high frequency dispersion is present. 

The plan of the paper is the following. In the next sec- 
tion we briefly review the equations governing the (1+1)- 
dimensional BEC in both the hydrodynamic limit and 
the dispersive case. In Sec. Ill we recall the brick wall 
model, and we show how one can use it in momentum 
space. In Sec. IV we apply the method to compute the 
entanglement entropy in the acoustic black hole with dis- 
persion. In Sec. V we discuss our results, and compare 
them with some numerical calculations. 



II. THE SETUP 

In the dilute gas approximation [l9|, the BEC can be 
described by an operator W that obeys the equation 



Bogolubov-de Gennes equation, we expand the field op- 



and we 



crator as <fi(t,x) = J^ - a,j(f)j(t,x 

find that the modes x) and (pj(t,x) satisfy the cou- 
pled differential equations [l4[ 
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where £ = K/ (mc) is the healing length of the condensate. 

In these settings, the dispersive effects are signaled by 
the presence of £, which depends on the local vel ocity 
of sound and that is a non-perturbative parameter [Tj] . 
Therefore, to study the case when dispersion is negligi- 
ble one must switch to the density-phase representation 
consisting in defining the density h\ and phase operators 



i\ via 
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along the lines of [12j . With these definitions, the limit 
£ — > is well defined and one finds the single equation 
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(6) 



thd t v = [-^ 2 + Kxt + g&$) $ , (i) 

where m is the mass of the atoms, g is the non-linear 
atom-atom interaction constant, and V cx t is the ex- 
ternal trapping potential. The wave operator satisfies 
the canonical commutation relations [^(t, x), ^(t, x')] = 
S 3 (x — x'). To study linear fluctuations, one substitutes 
W with Vo(l + </>) so that v satisfies the Gross-Pitaevski 
equation 



ihdt^a = ( -— V 2 + y cxt + gn ) v , (2) 
2m 



and the fluctuation cf> is governed by the Bogolubov-de 
Gennes equation 

i^=— V 2 +2^vU + m^+f), (3) 
2m \ W / 



where c = y^gn/m is the speed of sound, and n = 
I Wo 1 2 is the number density. We now focus on the 
(1 + 1)— dimensional case and consider a configuration 
with constant v and n, while the speed of sound c(x) 
smoothly decreases from the subsonic region to the su- 
personic one, and it is equal to v at x = 0, as in [l2j |. 

This is possible provided one modulates 5, and hence 
the speed of sound c, by keeping the combination gn + 
V ex t unchanged [l(|. In this way, Eq. ((2]) admits the 
plane- wave solution \&o = y/n e~Kp(ik[)X — iuj t) where 
V = hko/m is the condensate velocity. To study the 



The analogy with gravitational black holes comes about 
when one notices that the above equation can be written 
as 09 1 — 0, where the d'Alambcrtian is defined on the 
so-called acoustic metric 
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and where one assumes that 9\ does not depend on the 
transverse coordinates y and z [HI . The metric shows an 
event horizon located where c(x) = v and its structure 
is the same of the Painleve-Gullstrand line element, up 
to the conformal factor n/(mc). Note that Eq. © is not 
conformally invariant unless c is constant. 



III. THE BRICK WALL IN MOMENTUM 
SPACE 

We now turn the the calculation of the entropy with 
the brick wall method, originally described by 't Hooft 
in where it was applied to a (3+l)-dimensional 
Schwarzschild black hole. The method is based on the 
counting of the modes of a massive scalar field, with 
Dirichlet boundary conditions, defined inside a box of 
size L and placed at a distance e from the horizon. The 
result is that the entropy is proportional to the area of 
the horizon. Also, the proper distance between the hori- 
zon and the edge of the box is a constant of the order 
of the Planck length. The physical interpretation of this 
entropy is discussed in great detail in [||. 
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For (l+l)-dimensional gravitational black holes, the 
result is radically different and, at the leading order, the 
entropy reads 



(8) 



in the limit L/e » 1 [20(. In the case of the acoustic 
black hole studied here, the result is the same, at least at 
the leading order, if we do not account for high frequency 
dispersion. To show this, it is sufficient to assume that 
the modes that are solutions to Eq. © vanish at the 
boundaries of the segment [e,L], where e is located near 
the horizon, at x = 0, and L is the length of the near- 
horizon region, namely the region where we can linearize 
the speed of sound as c(x) ~ c(0) + kx on each side of 
the horizon. In this expansion 
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is the analog of the surface gravity of the black hole. The 
stationary solutions f(x) of Eq. ((6]) can be found with a 
WKB approximation by writing the solution as 



cxp 



f(x')dx' 



(10) 



This expression can be seen as the continuum limit of 
a function of the form exp(z^ n /c„) for n = 0...n max , 
where the sum counts the number of modes. Thus, in 
the continuum limit, the number of modes populating 
the interval [e, L] with frequency to is given by [5[ 

n(to) = \( f(x)dx~— In(-) . (11) 
irn J e ttk \ e J 

We now recall that the free energy and the entropy asso- 
ciated to massless spin-0 particles are given respectively 
by 



F = — 
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where E = hw and f3 is the inverse of the temperature of 
the black hole, p = {k B T)- 1 = 2-k/{Hk). 

By replacing the expression (TlUj) into Eq. © , and keep- 
ing only the leading terms in h we find that f(x) ~ 
Huj/(c(x) ± v). By expanding around the horizon, where 
c(x) ~ v + kx, we see that the dominating solution for 
small x is f(x) ~ Huj/kx, thus we find Eq. ©. This result 
confirms the validity of the brick wall model to compute 
the entanglement entropy when the acoustic metric has 
a horizon. In fact, the presence of a horizon is crucial to 
this result, as f(x) ocl/x precisely because of the vanish- 
ing of the gu term of the metric ([7]) at the horizon. We 
note further that Eq. (|8) agrees also with the calculation 
of the leading term of entanglement entropy in bi-partite 
spin-chains [2l|, including the non-critical case, where 
conformal symmetry is only approximate [22l |. as in the 
case studied here. 



To tackle the dispersive case, it is convenient to derive 
this result also in momentum space. One defines the 
Fourier transform 9\ of the modes via 



0i(x) 



dp 
v/2T 



l(p) 



so that Eq. ([5]) in momentum space becomes 
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where, in the near-horizon approximation, 
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The solutions f(p) to Eq. (fT4| can be written in terms 
of Whittaker function. However, for our purposes it is 
sufficient to use again the WKB method by substituting 
in Eq. (fl"4| the expression 
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: CXp 



Kp'W 



(16) 



If we consider the large vp/K limit, i.e. we select 
wavelengths much smaller than the near-horizon region 
(whose size is approximately v/k), and we recall that we 
are in the regime of linear dispersion uj — c(x)p , we find 
that, at the lowest order in the WKB expansion, 



tlUJ 
Hp 



(17) 



As we are in momentum space, we now count the modes 
with an associated momentum between p m i n and p ma x- 
The first value correspond the the largest wavelength ad- 
mitted in the near-horizon region and corresponds to the 
infrared contribution to the integral (jlip . The value p max 
is interpreted as the minimal distance that we can probe 
with our modes. As we are considering the hydrodynamic 
limit, this implies p m ax£ *C 1. In analogy with Eq. (jlip . 
the number of modes is defined as 



1 



f{p)dp 



V j I Pme 
TTK 



By following the same steps as above, we find that 



S=\\n 
6 



Pi, 



Pi, 



(18) 



(19) 



This expression is equivalent to Eq. (JSJ in terms of count- 
ing the numebr of degrees of freedom, however it stresses 
the fact that the entropy diverges in the ultraviolet non- 
locally, i.e. in no particular point in the near-horizon re- 
gion. In this respect, the expression above reflects more 
closely the properties of the entanglement entropy Q. 

Before considering the dispersive case, we recall that 
the near-horizon region has an extension L roughly given 
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by the speed of sound multiplied by the typical time At 
taken by a mode to cross this region. Therefore, by using 
Eq. (jHJ, we can write S — 1/(6 At). On the other hand, 
the surface gravity can be interpreted as the inverse of the 
time taken by a mode to cross the near-horizon reg ion. 
Therefore, we find that S = k/6, in line with |23(, up 
to a factor of two, that depends on considering bosons 
rather than fermions. Although there is no explicit time- 
dependence in the model, the Hawking mechanism still 
produces a steady flow of outgoing phonons at a rate 
that depends on k. Therefore, the time derivative of 
S should be interpreted as the production rate of the 
entropy associated to this flow. 



case. With the help of Eqs. (|12[) . we calculate the entropy 
in the form S = Si ea d + S COII where 

4 6 7H 6 \2nvpminJ 

Scon- = TTT^TT ° 2 ' (^) 

960 v z 

In these expressions 7 is the Euler constant and £(1, 2) is 
the first derivative of the ^-function evaluated at 2. To 
obtain this result, we assumed that the inverse temper- 
ature of the Hawking radiation is = 2tt/(Hk), i.e. it is 
not affected by dispersion [l5| . 



IV. DISPERSIVE CASE 



V. DISCUSSION 



We now consider the dispersive case, and evaluate 
the contribution to the entropy given by high frequency 
modes. The two equations of the system Q can be easily 
decoupled in momentum space. By defining \E' ± (i,a;) = 
cxp(iujt)[4>(x)iip(x)], we find that the Fourier transforms 
of Eqs. Q reads 



c 2 *+(p) - 



(w — pv) h 2 p 



4m 2 



*"(p) 



2m{u-pv) ^ 

hp* 



. (20) 



The first of these equations can be solved with the WKB 
method in the near-horizon approximation c = v + ind p . 
At the leading order, we find 



(21) 



where £0 = £(- T = 0) is the healing length at the horizon. 
The number of modes is 




h(ui) 



dx 



x > 



x" 

T 



, (22) 



where x = p^o and a = lj^q/v. The integration bound 
aries are fixed by the positivity of p and of the argument 
of the square root, namely < x < 
1 



1 + \/T+ 2a U 
However, as the wave- 



VI - 2a < x < 1 + VI - 2a 
lengths of the physically realistic modes, given by v/uj 
must be much larger than £oj we see that a« 1. Thus 
we can simplify the above integral by expanding the in- 
tegrand function and the upper integration limit as 



h{fjS) 



tt£ k 



dx 



x 
8^ 



(23) 



where e = ^oPmin C a is set to cope with the same 
logarithmic divergence encountered in the non-dispersive 



We first note that the correction S^orr to the entropy 
is set by the square of the size of the near-horizon region 
L ~ v/k measured in units of the healing length. Re- 
markably, this term is very similar to the one found in 
the case of the one dimensional spin- 1/2 Hciscnberg XX 
chain in a magnetic field [24[ . 

Another important observation is that the leading term 
S\ c &d is no longer completely arbitrary as in the hydro- 
dynamic case. In fact, dispersion affects the integration 
boundaries in Eq. (|22[) . which are no longer put by hands, 
but are fixed in the ultraviolet by the system itself, as the 
dispersion relation in this system is fixed by the proper- 
ties of the BEC and can be found to be expressed by the 
function UM 



(lj - vp) 2 = c 2 [p 2 + 



e P 4 



(26) 



As a result, the leading term turns out to be a numeri- 
cal constant plus a logarithmic correction that depends 
on the infrared regulator. Physically, the argument of 
the logarithm k/ (2irvp m i n ) can be written as the ra- 
tio of the maximum wavelength allowed in the system 
Amax = 1/pmin and the size of the box, i.e. of the near- 
horizon region, L. It is reasonable to expect that, by 
taking in account the dependence of 0\ upon transverse 
directions (while keeping the potential a function of x 
only) one can provide for an effective mass able to act as 
an infrared regulator. This issue will be investigated in 
a future paper. 

To further check our results, we compute the entropy 
numerically, and we cope with the infrared divergence in 
two ways. In the first, we set the value of x m in = 10~ 12 
in the integral (f2"2")l . In the second, we subtract to the in- 
tegrand the function a/x and we let x m i n = 0. The upper 
integration limit is set at x max = 10 -4 in both cases. The 
two results are plotted in Fig. ([T]), and we see that the en- 
tropy is in fact constant in terms of the normalized tem- 
perature £fikB,T/{vh) ~ £o/(27rL). The different numeric 
values of the plateaux depend only on the choice of the 
infrared regularization. If one considers typical experi- 
mental values for Rubidium atoms, like v = 4x 10~ 3 m/s, 
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£o = 2 x 10~ 7 m, k = 2, 7 x 10 3 HzAhe Hawking temper- 
ature is of the order of few nK [12|], which corresponds 
to values on the horizontal axis around 0.05. We see 
that the entropy is constant in a large range containing 
this value. This is expected as the first order correction, 
which can be written as S'corr — — (7r£o/cBT/4i>fi,) 2 /15, is 
very small for these values. 

In summary, we have verified that the scaling behavior 
of the entropy in the hydrodynamic limit is the same as 
the one predicted by conformal field theory, by using a 
method inspired by the brick wall model for astrophysi- 
cal black holes. When dispersion is taken into account, 
we found a correction that is similar to the one calcu- 
lated for the entanglement entropy of certain spin-chain 
systems. Also, the leading term appear to be a constant 
determined uniquely by the infrared cut-off. These ele- 
ments, although not a rigorous proof, strongly support 
the interpretation of the brick wall entropy as due to the 
entanglement of the phonon pairs created via the Hawk- 
ing mechanism. 
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FIG. 1: Normalized entropy versus the temperature for small 
temperatures. The bottom curve is obtained by subtracting 
the function a/x from the integrand of Eq. I|22|l and setting 
Xxnin = 0. The top curve is obtained by setting x mln = 10~ 12 . 
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